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An extension of the Gutzwiller trace formula is given that includes diffraction effects due to hard 
wall scatterers or other singularities. The new trace formula involves periodic orbits which have arcs 
on the surface of singularity and which correspond to creeping waves. A new family of resonances 
in the two disk scattering system can be well described which is completely missing if only the 
traditional periodic orbits are used. 



Gutzwiller's trace formula ||l| is an increasingly popu- 
lar tool for analyzing semiclassical behaviour. Recently, 
it has been demonstrated that using proper mathemat- 
ical apparatus, like Gutzwiller- Voros zeta functions, 
cycle expansions 1^ or Quantum Fredholm Determinants 
[^, the trace formula can successfully predict individual 
resonances of open scattering systems Jp], it is possible 
to compute systematic h corrections p] and to extend 
the formula for systems with regular domains |7J. The 
physical content of the trace formula is the geometrical 
optical approximation of quantum mechanics via canon- 
ical invariants of closed classical orbits. This approxi- 
mation is very accurate when periodic orbits sufficently 
cover the phase space of the chaotic system. This is not 
the case when the number of obstacles is small or their 
distance is large compared to their typical size. In such 
cases it is very important to take into account the next- 
to-geometrical effects. 

In this paper we study how the Geometric Theory of 
Diffraction (GTD) for hard core potentials can be incor- 
porated in the periodic orbit theory. Such a problem oc- 
curs where the wave length of a quantum mechanical (or 
optical) wave is very large compared with the spatial vari- 
ation of a repulsive potential, e.g. at the boundaries of mi- 
crowave guides, optical fibers, superconducting squids, or 
circuits in the ballistic evolution of electrons, i.e. in most 
of the devices used for so-called macroscopic quantum 
mechanical (or optical) experiments. First we summa- 
rize from the classical papers of J. B. Keller ||l^ how the 
Green's function in the shadowed regions of configuration 
space can be computed by the concept of GTD. Then we 
incorporate the periodic rays with diffracted ray arcs to 
the trace formula. It is possible to derive the diffraction 
part of the Green's function directly from the semiclas- 



sical approximation of the Feynman path integral in the 
neighborhood of a hard curved wall [g| . We have choosen 
Keller's approach here, since it is more suitable for com- 
puting the trace of the Green's function, without further 
approximations. 

It has been known for quite a long time that the scat- 
tering amplitude in the shadowed region behind an ob- 
stacle can be well reproduced by allowing diffracted rays 
in addition to geometrical ones |Q . The diffracted rays 
connecting two points in the configuration space in the 
presence of sharp objects can be derived from an ex- 
tension of Fermat's variational principle of classical me- 
chanics [|0|. The usual Fermat principle states that the 
classical trajectories connecting two positions <7_4 and gg 
in configuration space are those smooth curves which 
make the action stationary. If the configuration space 
is bounded by hard walls, a generalized variational prin- 
ciple, introduced by J. B. Keller [nO|, has to be applied. 
This principle requires new classes of curves: We have to 
consider for each triplet of integers r, s, t > the class of 
curves Vrst with r smooth arcs on the surface, s points 
on the edges and t points on the vertices of the bound- 
ary or the discontinuity. The curves of the GTD are 
those which make the action stationary within one of the 
classes Vrst ■ The class I?ooo corresponds to the usual ge- 
ometrical orbits. In this paper, we concentrate on two di- 
mensional problems, where the simplest nontrivial curves 
are of class I?ioo and I?ooi (Fig- la and b respectively), 
whereas edge diffraction I?oio is not possible. Once we 
know the generalized ray we can compute semiclassically 
the Green's function G((7^, gg, -E) tracing the ray along 
it ||l^. In Fig. la (I'loo) the trajectory - obtained from 
the generalized variational principle - is tangent to the 
surface of the hard wall obstacle at the points A' and 



B'. The Green's function G{q^, q^',E) can be computed 
semiclassically by the energy domain Van-Vleck propa- 
gator (for a single classical trajectory) 



G{q,q\E) = 



27r 



i27rih)y^ 



i^y^(g,g',i?)e^^('-''^)- 



(1) 



where Dv{q,q',E) = \det{-d^S/dq,dq'j)\/\q\\q'\ is the 
Van-Vleck determinant and v is the Maslov index (see 
ref. |1^ for definitions). 

When the geometrical ray hits the surface of the obsta- 
cle it creates a source for the diffracted (creeping) wave. 
The strength of the source is proportional to the Green's 
function at the incidence of the ray 



^diff 



DG. 



(2) 



The diffraction constant D depends on the local geometry 
and the nature of the diffraction. It has been determined 
in Ref [ |lO| from the asymptotic semiclassical expansion 
of the exact solution in some simple geometry [Q (see 
also fO). Its form is 



D, = 2I/33-2/3 






(3) 



Here Ai'{x) is the derivative of the Airy integral Ai{x) = 
Jp dt cos{xt—t^), k = yj2mE/h is the wave number, p is 
the radius of the obstacle at the source of the creeping ray 
and xi are the zeroes of the Airy integral, which can be 
approximated as xi = Q^^^ (37r(/— l/4))^/'^/2 in the semi- 
classical limit. The index I > 1 refers to the possibility of 
initiating creeping rays with different modes, each with 
its own profile. In practice only the low modes contribute 
to the Green's function. The source then initiates a ray 
creeping along the surface. During the creeping of the 
ray the amplitude decreases, which can be understood as 
a process analogous to the radiation processes of electro- 
dynamics. The radiated intensity is proportional to the 
intensity of the ray: 



—Ai{s,Ef=-2ai{s,E)Ai{s,E)' 
as 



(4) 



where s is the length measured along the surface and 
Ai{s,E) is the complex amplitude of the Green's func- 
tion along the surface. The coefficient ai{s,E) depends 
on the local curvature of the surface, 1/ p{s), and it has 
the structure ai{s,E) = a;;e~"/'^(A;/6p(s)2)i/3 (gee Ref. 
|l3| ) , where the index I refers again to the different modes 
of the creeping wave. The Green's function for the creep- 
ing ray of mode I is then given by 



Gf{qA'.qB',E) 



-j^dso,l(s,E)l.s(q^,,q^,,E) 



(5) 



where L is the length of the arc of the creeping ray, and 
S{qA',qB',E) is the action along it. The creeping ray 



at the point B' initiates a pure geometrical ray. The 
source of this ray is located in B', and its strength is 
again given by equations (||) and (||) due to the invari- 
ance of the Green's function against the interchange of 
the variables qA> and qs'- The total Green's function is 
then the product of the Green's functions and diffraction 
coefficients along the ray: 

oc 
G{qA,qB,E)^G{qA,qA':E)Y,Di,A'Gf{qA',qB',E) 



1=1 



xDi^B'G{qi3',qB,E). 



(6) 



In a general situation, when the ray consists of several 
pure geometric and creeping arcs, the Green's function 
can also be written as a product of partial Green's func- 
tions and diffraction constants. 

To incorporate diffraction effects into the trace for- 
mula, one should compute the trace of the Green's func- 
tion derived above. As in the case of the Gutzwiller trace 
formula - derived from a pure geometrical approxima- 
tion of the Green's function - the trace receives the lead- 
ing contributions from tubes encircling the closed curves, 
which now can have diffractional arcs too. We can handle 
separately the pure geometric cycles and the cycles with 
at least one diffraction arc along one of the obstacles: 



Tr G{E) w Tr Gg{E) + Tr Gd{E), 



(7) 



where TtGg{E) is the ordinary Gutzwiller trace formula, 
while TrGoiE) is the new trace formula corresponding to 
the non-trivial cycles of the GTD. TtGdJE) can be com- 
puted by using appropriate Watson |l4yil| ] contour inte- 
grals. For technical details we invite the reader to Refs. 
|l5| and iQ. Here we communicate the general result 
and the detailed calculation will be published elsewhere 
]l7| . If we denote by g^, i = 1, . . . , n (with qn+i = qi) the 
points along the closed cycle, where the ray changes from 
diffraction to pure geometric evolution or vice versa (see 
A' and B' in Fig. la), the trace for cycles with at least 
one diffractional arc can be expressed as the product 



TrGniE)^ J2 '^l[D{q,)G{q„q,+i,E), 



(8) 



Cycles 



where T(E) is the time period of the cycle (without re- 
peats) and D{qi) is the diffraction constant (^) with the 
radius of curvature p given locally at the point qi. The 
creeping mode index I and the corresponding summations 
(see e.g. (0)) are surpressed here for keeping the nota- 
tion simple. G{qi, qt+i, E) is alternatingly the Van-Vleck 
propagator, if qi and g^+i are connected by pure geo- 
metric arcs, or is given by dq) in case qi and qi+i arc the 
boundary points of a creeping arc. Note this formula does 
only apply for cycles with at least one creeping section. 
Such cycles have the special property that their perti- 
nent energy domain Green's functions are multiplicative 



(see e.g. (^ with the summations over the creeping mode 
numbers of course included consistently) . This does not 
hold for pure geometrical cycles. 

The eigenenergies can be recovered from the 
Gutzwiller-Voros spectral determinant O] A{E), which 
is related to the trace formula as 



TTG{E)^—\nA{E). 



(9) 



The full semiclassical determinant can be written as the 
formal product of two spectral determinants, one cor- 
responding to pure geometrical and one to new cycles 
A(£') = Ag(^)Ai3(^) due to the additivity of the 
traces. The product is only formal, since the eigenen- 
ergies are not given by the zeros of Ag{E) or Ad{E) 
individually, but have to be calculated from a curvature 
expansion of the comfomed determinant A{E) itself. 
The diffraction part of the spectral determinant is 

An{E)^e^p(-Y, -Y[[D{q!)G{qf,q!^,,E)Y), 



r 
=1 i=i 



(10) 

where the summation goes over closed primitive (non- 
repeating) cycles p and the repetition number r. The 
product of Green's functions should be evaluated for q^ 
belonging to the primitive cycle p. After summation over 
r, the spectral determinant can be written as 



with 



AD{E)=\{{l-tp) 



\{D{ql)G{q^,ql^,,El 

4=1 



(11) 



(12) 



where gf belongs to the primitive cycle p. As mentioned 
before the mode numbers I of the diffraction constants 
and the corresponding summations have been surpressed 
for notational simplicity; they can be easily restored as 
e.g. in the final expression ([l^). 

We can conclude that the diffractional part Ad(£') of 
the spectral determinant shares some nice features of the 
periodic orbit expansion of the dynamical zeta functions 
PI, and it can be expanded as 



AD{E) = l-Y,h + ^tph 



(13) 



p,p' 



The weight (n2) has the following property which helps in 
radically reducing the number of relevant contributions 
in the expansion: If two different cycles p and p' have at 
least one common piece in their diffraction arcs, then the 
two cycles can be composed to one longer cycle p + p' 
and the weight corresponding to this longer cycle is the 
product of the weights of the short cycles 



Ep+p' 



^p ' '^p' • 



(14) 



As a consequence, the product of primitive cycles, which 
have at least one common piece in their diffraction arcs, 
can be reduced in such a way that the composite cycles 
are exactly cancelled in the curvature expansion 



l[{i-tp) = i-Y,tb, 



(15) 



where tt are basic primitive orbits which can not be com- 
posed from shorter primitive orbits. 

To demonstrate the importance of the diffraction ef- 
fects to the spectra, we have calculated the Ai resonances 
of the scattering system of two equally sized hard circular 
disks with disk separation R — 6a, where a is the radius 
of one disk. In this system there is only one geometri- 
cal periodic cycle along the line connecting the centers 
of the disks. Its stability Ap = 9.8989794 and action 
Sp = kLp = k ■ 4a yield the geometrical part of the spec- 
tral determinant |18| , [16[ 



AG{k)^llil 



ikLj} 



A 



(l+4i)/2 



(16) 



where k — \/2mE/h and 2m ^ h — 1, and leads to the 
following predictions for the semiclassical Ai resonances 



fc;<;5 = (7r(2n ~ 1) - i^ lnAp)/Lp 



(17) 



with n — 1,2,3,---. Note in the above expressions 
(1 + 4j)/2 replaces the usual weight (1 -I- 2j)/2, since 
the geometrical orbit in the two-disk problem lies on the 
boundary of the fundamental domain p8| . 

Fig. Ic shows the first four new basic cycles in the 
fundamental domain p8| . We computed the geometrical 
data of the first ten orbits and used them to construct the 
creeping and geometrical Green's functions. The semi- 
classical Green's function in free space is asymptotically 
(/ci?>l) 



Go{q,q',E) = - 



4 \nkR 



1/2 



ikR~i-r 

e 4, 



(18) 



where R = \q — q'\. If the ray connecting q and q' is 
reflected once or more from the curved hard walls before 
hitting tangentially one of the surfaces, we can keep track 
of the change in the amplitude by the help of the Sinai- 
Bunimovich curvatures. 

By computing the curvatures Ki right after the reflec- 
tions, and knowing the distances U between the z-th and 
the [i + l)-th points of reflections, the factor R in the 
Green's function (|8|) has to be changed to the effective 
radius R^^ — Rq YlTLii^ + hi^i) where Rq is the distance 
between q and the first point of reflection along the ray 
starting from q, and m is the number of reflections from 



a disk. The effective radius R-l^, the length of the geo- 
metrical arc Lf and the length of the diffraction part L^ 
of the first ten orbits with creeping sections are listed in 
Table |. To each cycle in the list, there is a whole se- 
quence of cycles which wind around the disk TOw times. 
For these orbits one has to add 2Tram^ to the diffraction 
length Lj^. The diffraction part of the spectral determi- 
nant is finally given by 



Ao(fc) = l-^(-inC, 



b,l 



^i/ogi/./ ±^gi>>,vijj, -r-i^b 



k^/^xlRf 



1 



1 __ p27r{ik — ai)a ' 



(19) 



where Ci = 7r^/^3^'*/^2^^/^/yli'(a;/)^, and the summation 
for the windings TOw gives the factor 1/(1 — e^'^(*'^~"')°). 
We computed the spectra by truncating the product 
Ac{k)AD{k) at maximal cycle length 5 and using only 
the I = 1 term in the now restored summation over the 
creeping mode number. The exact quantum mechanical 
resonances were computed following Ref. |15|] . 

The leading semiclassical resonances are given equally 
well with and without creeping modifications. In Fig. 2 
we can see that the new formula describes the resonances 
of the two disk system with a few-percent error, while 
the computation based on the geometrical cycle alone, 
Eq. (uw, gives completely false results for the next-to- 
leading resonances (see Eq.(n^)). 

The authors are very grateful to P. Cvitanovic for en- 
couraging the project. A.W. thanks NORDITA and NBI 
for hospitality during his stay. G.V. and P.E.R. thank 
SNF for support and G.V. acknowledges the support of 
the Szechenyi Foundation, Phare Accord H 9112-0378 
and OTKA F4286. 



[10] J. B. Keller, J. Opt. Soc. Amer. 52 116 (1962) J. B. Keller, 
in Calculus of Variations and its Application (American 
Mathematical Society, 1958) p. 27 B. R. Levy and J. B. 
Keller, Comm. Pur. Appl. Math. XII, 159 (1959) B. R. 
Levy and J. B. Keller, Cann. J. Phys. 38, 128 (1960) 

[11] W. Franz, Theorie der Beugung Elektromagnctischer 
Wellen, Springer Verlag, Berlin (1957); Z. Naturforschung 
9a, 705 (1954) 

[12] M. C. Gutzwiller, Chaos in classical and quantum mechan- 
ics, Springer Verlag, New York (1990) 

[13] B. Schrempp, F. Schrempp, Nucl.Phys.B163,397 (1980) 

[14] G. N. Watson, Proc. R. Soc. London A95, 83 (1918) 

[15] A. Wirzba, CHAOS 2, 77, (1992) 

[16] A. Wirzba, Nucl. Phys. A560, 136, (1993) 

[17] G. Vattay, A. Wirzba and P. E. Rosenqvist, in prep. 

[18] P. Cvitanovic, B. Eckhardt, Nonlinearity 6, 277, (1993) 



TABLE L Table of the first ten basic cycles tt which in- 
clude creeping sections in the fundamental region of two-disk 
problem (with disk separation 7? = 6a). The cycles are la- 
belled by their number nif, of geometrical reflections from one 
of the disks. The length of the geometrical arc Lf, the effec- 
tive radius Rl and the length of the diffraction section L^ 
are listed in units of the disk radius a. 



nib 


L°/a 


Rf/a 


L?/a 





5.656854249492 


5.656854249424 


3.821266472498 





6.000000000000 


6.000000000000 


3.141592653589 


1 


9.832159566199 


58.16784043380 


3.476488812029 


1 


9.797958971132 


58.78775382679 


3.544308495170 


2 


13.81654759452 


578.1406653460 
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17.81499162871 


5729.649817456 


3.510488616089 
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17.81464272590 


5732.235502463 


3.511180541615 


4 


21.81483475355 


56728.70010470 


3.510799703655 


4 


21.81479950722 


56732.26871144 


3.510869602322 



FIG. 1. The simplest classes Oioo (a) and ©ooi (b) of curves 
in two dimensions. In the window (c): the first four basic 
orbits in the fundamental domain of the two-disk system. 
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FIG. 2. Resonances for the Al subspace of the two-disk 
system (with disk separation R — 6a) in the complex k plane 
in units of the disk radius a. The diamonds label the exact 
quantum mechanical resonances, which are the poles of the 
scattering matrix. The crosses are their semiclassical approx- 
imations including the diffraction terms derived in this paper. 
The boxes refer to the ordinary Gutzwiller semiclassical ap- 
proximation (eq.(17), j = 0,1), where the diffraction effects 
are not included. 
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